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ABSTRACT 
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(Nj , Recent Spitzer infrared measurements of hot Jupiter eclipses suggest that 

eclipse mapping techniques could be used to spatially resolve the day-side pho- 
^ ! tospheric emission of these planets using partial occultations. As a first step in 
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this direction, we simulate ingress/egress lightcurves for three bright eclipsing hot 
CN \ Jupiters and evaluate the degree to which parameterized photospheric emission 

\ models can be distinguished from each other with repeated, noisy eclipse mea- 

surements. We find that the photometric accuracy of Spitzer is insufficient to 
use this tool effectively. On the other hand, the level of photospheric details that 
could be probed with a few JWST eclipse measurements could greatly inform 
hot Jupiter atmospheric modeling efforts. A JWST program focused on non- 
parametric eclipse map inversions for hot Jupiters should be actively considered. 

2 ■ 1. Introduction 

CO 

Eclipses can be used as powerful tools to spatially resolve the photospheric emission 
properties of astronomical objects. During ingress and egress, the partial occultation effec- 
tively maps the photospheric emission region of the object being eclipsed. With arbitrarily 
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high photon statistics, the (spatial) resolving power attainable with such an astronomical 
telescope is unrivaled. In the past, eclipse mapping methods have th us been used with g reat 



success to measure limb darkening laws on stellar companions (e.g. iWarner et al.lll97ll ). to 
map the surface appearance of Pluto (see review by lSternlll992l ). or to confirm the Keplerian 
nature of accretion dis ks in Cataclysmic Variables, on the basis of their radial emission pro- 



files (e.g. lRuttenlll996l). Given the recen t Spitzer measurements of infrared eclipses for three 



hot Jupiters (jDeming et al. 2005 . 2006 ; Charbonneau et al. 20051). it i s natural to expL 



the potential for eclipse mapping of hot Jupiters ([Williams et al 



lore 



20061 ). 



The regime of atmospheric circulation present on hot Jupiters may be unlike any of 
the familiar cases encountered in the Solar System. Hot Jupiters are massive, gaseous giant 
planets found in close, circul ar orbits around th eir parent stars, with orbital periods on 

2006h. General arguments suggest th at these 



the order of a few days 



planets are tidally locked ( Rasio et al. 19961 ; 



e^ 



Butler et al 



Lubow et al.lll997l : lOeilvie & Linll2004h . so that 



their permanent day-sides are continuously subjected to intense stellar irradiation, while 
their night-sides may be heated only by a much more modest internal energy flux. In the 
presence of such an uneven energetic forcing, leading to significant (horizontal) temperature 
gradients, atmospheric winds will attempt to redistribute the heat more evenly around the 
planet. The efficiency of this horizontal heat redistribution is an important open question 
for hot Jupiters as it largely determines their observational properties. Several attempts 
have been made to address this general circulation problem for hot Jupiter atmospheres 



using different approache s (jShowman &; Guillotl l2002t ICho et al.l 12003 ; iMenou et al.l 12003 



Cooper Showman! 120051 . see also Burkert et al. 2005), and interesting constraints are now 
being placed on these models through the recent infrared measurements. 

Three hot Jupiters have been directly detected in i nfrared with the S pitzer Space 
Te lescope: HD189733b (fPeming et al-l boOffl. HD209458b JDeming et al.ll2005h . and TrES- 
1 (ICharbonneau et al.l 120051 ) . The planetary day-side flux is deduced from the secondary 
eclipse depth, when the planet is hidden behind its star. Based on such flux measurements 
at different wavelengths, several groups have computed spectral emission models for various 
prescribed levels of heat redist ribution and concl uded that at least some amount of heat 



redistribution must be present (ISeager et al 



Barman et al. 2005; Burrows et al. 



2006 



2005 



Burrows et al.ll2005l : iFortnev et all 12005 



Fortney et al.ll2006l ). These secondary eclipse con- 



straints remain weak at the moment but they could rapidly improve in quality in the future. 
In addition, a first hot Jupiter infrared orbital phase curve for t he non-eclipsing planet Ups 
And b has been reported very recently ([Harrington et al.l 120061 ) . The large inferred ampli- 
tude of the phase curve, measuring the day- to night-side temperature difference, suggests a 
rather weak role for heat redistribution in this case. More work is likely needed to interpret 
the phase curve reliably, and in particular to address the complications that may arise from 
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the e xistence of different regimes of heat redistribution at different heights in the atmosphere 



e.g. 



Seager et al.ll2005l : llro et al.ll2005l : ICooper fc Showman! 120051 ). 



This recent progress in direct infrared detections of hot Jupiters leads us to consider 
the possibility that eclip se mapping techniques could also be used in these systems (see 
also IWilliams et al.ll2006l ). Contrary to existing measurements, which only constrain disk- 



integrated quantities, eclipse mapping of a hot Jupiter can provide a "spatially resolved 
map" of its day-side photospheric disk, with "pixel" shapes and sizes determined by eclipse 
geometry and photometric accuracy. In principle, one can measure the emission of any 
individual "pixel" directly from the shape of the ingress and egress lightcurves, without 
making any significant assumption on the nature of the photospheric emission itself. Such 
"non - parametric" inversion techniques can be used in the eclipse mapp ing context (e.g. lHorne 
198,4 ISpruitlll994l : lRuttenlll996l . Il998l : iBohingej 1200(1 lBaptistall2~00lh . but in this first study 
we focus on a much simpler parametric method to evaluate the accuracy with which details of 
the planet's day-side photospheric emission can be constrained with Spitzer now and JWST 
in the future. Specifically, we use physically-motivated, parameterized models to describe the 
photospheric emission properties of hot Jupiters. We generate artificial eclipse curves, with 
appropriate levels of noise, for the three detected eclipsing hot Jupiters and we determine 
how well the various models for their day-side photospheric emission can be distinguished 
from each other with repeated eclipse measurements. This study is thus strongly model- 
dependent. As will be emphasized throughout this paper, it will be important in the future 
to perform non-parametric studies of eclipse mapping for hot Jupiters, since this would be 
a much more reliable way to interpret actual data. 



Williams et al.l (120061 ) first mentioned the possibility of resolving the "surfaces" of hot 
Jupiters with secondary eclipses. These authors simulated secondary eclipse lightcurves and 
introduced the concept of a "uniform time offset" as a useful observable to differentiate the 
lightcurves of various photospheric emission models from that expected for spatially uniform 
emission. In particular, they find that Spitzer should be able to distinguish betw een the 
emission properties predicted by the circulation model of ICooper fc Showman (120051) and the 



case o f spatially uniform emission. Our work shares similarities with that of IWilliams et al. 



( 120061 ) but it also differs from it in a number of ways, most significantly in that we use the 
full shapes of ingress/egress curves as observational diagnostics. 

The set of photospheric emission models considered in our study are described in Sec- 
tion [2j Section [3] explains our method to simulate eclipse curves and generate the associated 
observational noise, as well as the statistical tools we use to analyze the data set resulting 
from repeated eclipse simulations. Section 0] details the results of our parameter space survey 
for Spitzer, while Section [5] focuses on JWST. We summarize our main results and propose 
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extensions of this work in Section [61 



2. Photospheric Emission Models 

Throughout this analysis, we postulate the existence of a well-defined photospheric 
surface for any spectral band of observation, so that the problem of generating ingress/egress 
lightcurves is reduced to that of specifying temperatures and/or emission fluxes at various 
locations on the planet's photospheric disk. 

At the simplest level, there are two extremes for horizontal heat redistribution in hot 
Jupiter atmospheres: complete or none. In the case of complete heat redistribution, presum- 
ably achieved by atmospheric wind and wave transport, a nearly constant temperature across 
the planetary disk is to be expected (modulo limb darkening effects). In the absence of any 
heat redistribution, however, a strong temperature gradient (set by purely local radiative 
equilibrium) is expected on the day-side, while the night-side would remain uniformly very 
cold under the influence of the modest internal energy flux. In between these two extremes, 
photospheric disk emission properties would be set by a complex non-linear combination of 
radiative and advective processes in the atmosphere. Two groups have performed global nu- 
merical simulations exploring the circulation regime in hot Ju piter atmospheres. The multi- 
dimensional model presented by ICooper fc Showman! (120051 . see also Showman & Guillot 
2002) predicts that winds could advect the hottest part of the atmosphere ~60° in longi- 
tude downstream of t he sub- stel l ar po int. The turbulent shallow-layer models presented by 
Cho et al.l (120031 ) and lCho et al.l (120061 ) predict the emergence of temporally variable temper- 
ature structures associated with two large-scale polar vortices, revolving around the planet. 
The observability of such features, however, was found to depend on a combination of global 
wind strength and net radiative forcing in the atmosphere. In particular, a situation ap- 
proaching pure radiative equilibrium can be reached in t he limit of weak a tmospheric winds 
and transport. For a detailed parameter space study, see I Cho et al.l (120061 ) . 



These various scenarios provide us with four distinct and physically well-motivated 
models for hot Jupiter day-side photospheric emission, from which we can build synthetic 
ingress/egress lightcurves. We label and parameterize these four scenarios as follows: 



RadEq Model: The radiative equilibrium model corresponds to the no heat redistribution 
case, for which one expects a day-side temperature profile T cx cos 1//4 (#), where 9 is 
the angle away from the sub-stellar point. This arises simply from the geometry of 
stellar irradiation on the day-side. We assume a constant temperature night-side. The 
specific functional form that we adopt for the photospheric temperature is T(9) = 
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[T 4 cos(#) + T n 4 ight ] 1/4 for < 6 < 90° and T = T night for 6 > 90° (on the night- 
side). We calculate the radiative equilibrium temperature at the sub-stellar point as 
T = T* y/B^/a, where T* is the effective temperature of the host star, is the stellar 
radius, and a is the orbital semi-major axis of the planet. 

Uniform Model: In this model, complete, efficient heat redistribution is assumed, re- 
sulting in a constant temperature everywhere at the photosphere. We adopt T = 



4 _|_ rp4 

p ' night 



1/4 



Here the "redistributed" equilibrium temperature is calculated as 



T P = T*^RJ2^{=T /V2). 



CS05-like Model: This is a parameterization of the atmospheric model of lCooper fc Showman 
(j2005l ). in which a hot spot is shifted by ~60° in longitude downstream of the sub-stellar 
point. For the hot spot itself, we adopt the same functional form as in the RadEq model 
above, except that the hottest region is only ~400K hotter than the cold "night-side" 
(and shifted away from the sub-stellar point by 60° in longitude). The left panel in 
Figure [1] shows the corresponding temperature map partially eclipsed. 

Cho03-like M odel: This is a sn apsh ot tempera t ure m ap taken from one of the atmospheric 
models of Icho et all j^OoJ and lcho et all hood ) for HD209458b. Two cold spots as- 
sociated with circumpolar vortices rotate around the planet semi-periodically, with 
a period typically several times that of the planet itself. The specific shallow-layer 
model parameters in the simulation used here are: global wind speed U = 400 m s _1 
and relati ve amplitude of t hermal forcing rj = 0.2 (for a full discussion of the parameter 



space see iCho et al.ll2006l ). The cold spots are ~10% cooler than the average atmo- 
spheric temperature. A fixed phase of the atmospheric temperature pattern has been 
chosen here for simplicity (it rotates around the planet with time), approximately in 
such a way that the greatest difference in ingress/egress lightcurve shape is achieved 
relative to the RadEq model lightcurve shape. The right panel in Figure [1] shows the 
corresponding temperature map partially eclipsed. 



We expect these parameterized, physically-motivated models to cover a wide enough 
range of possibilities for us to be able to assess the level of photospheric details that 
can be probed via the corresponding ingress/egress lightcurves. In particular, in the or- 
der RadEq/CS05-like/Cho03-lik3/Uniform, these models constitute a sequence of increasing 
deviation from simple radiative equilibrium, towards photospheric temperature uniformity. 

We use these photospheric emission models to simulate ingress/egress curves for the 
three, bright hot Jupiters with Spitzer eclipse measurements: HD 189733b, HD209458b, and 
TrES-1. All the system parameters used in our analysis are listed in Table 1. The range of 
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temperatures in each photospheric model is adjusted to appropriate values for each planet 
of interest, using that system's parameters. From the above description of the RadEq and 
Uniform model, it is clear that we have assumed negligibly small albedos in all cases. We 
find tha t a change of Bond albedo from to 0.3 (a range characteristic of most Solar System 
planets; Ide Pater fc Lissauerll200ll ) has little influence on our main results, considering other 
systematic effects present. For computational reasons, we set the value of T nig ht to 700 K in 
all our models. Here again, we find that this choice is not crucial because very little of the 
night side is generally visible to a distant observer, given the tidally-locked eclipse geometries 
considered here. 



From the detailed model of ICooper fc Showman! (120051 ). values of T and T night for 
HD209458b's photosphere are approximated as 1450K and 1010K, respectively. In our 
CS05-like models for the two other hot Jupiter systems, these values are rescaled in lin- 
ear proportion to the value of the "redistributed" equilibrium temperature of the planet of 
interest ( that is T p = T^y / R*/2a for zero albedo; see Table 1). In the detailed HD209458b 
model of ICho et al.l (120031 . 120061 ) used here, the hottest and coldest temperatures reached 
overall are approximately 1460K and 1230K, with an average temperature of 1390K. In our 
Cho03-like models for the two other hot Jupiter systems, these values are also rescaled in 
linear proportion to the value of the "redistributed" equilibrium temperature of the planet 
of interest. 

The conversion from local photospheric temperature to emission contributed by each 
apparent area element of the planetary disk (see §3. II for projection details) is performed in 
two different ways. In a set of blackbody emission models, a simple Planck function is used. 
In a second set of models, a grid of detailed atmospheric spectral models derived from plane- 
parallel calculations is used to evaluate the local emission for effective temperatures in the 
range 700-2000K. No temperatures outside of this range of validity for the spectral models 
were allowed in our photospheric emission models. The plane-parallel spectr al models are 
essentially extensions of the clo udless models describ ed in ISeager et al.l (120051 ) . Contrary to 



more detailed calculations (e.g. iBarman et al.ll2005l ). we ignore the 2D geometry of the ra- 



diative transfer problem when we use these plane-parallel calculations and implicitly assume 
isotropic emission of each planetary surface element in our photospheric emission models. 

Planetary limb darkening (or brightening) could in principle be important in hot Jupiter 
atmospheres. It is not included in our basic models given the very simple method adopted 
(which assumes a well-defined photospheric temperature field and an isotropic emission). 
To estimate the effect that strong limb darkening could have on our overall analysis, we 
have included an optional limb darkening function in our emission models which can be 
applied to any of the models under consideration. The functional form adopted for this 
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strong limb darkening is chosen to correspond to a geometrical decrease of the emitted 
flux with the cosine of the angle away from the sub-observer point, (i.e. F = F o cos(0), 
where <ft is the angle away from the sub-observer point). Since our analysis relies on well- 
defined photospheric temperatures, we have implemented this in our "limb-darkened" models 
with a simple decrease of the photospheric temperature as cos 1 / 4 (0), superimposed on the 
photospheric temperature field of any of the four emission models defined above. 

Another systematic source of errors for hot Jupiter eclipse mapping studies could be 
related to finite eclipse timing uncertainties. Even for planets with measured orbital eccen- 
tricities consistent with zero, the presence of an as yet undetected pertu rbing companion 



could lead to significant timing variations from one eclipse to the next (lAgol et al.l 12005 



r 



Holman fc Murray! 120051 ) . Careful measurements of transit times before and after any sec- 
ondary eclipse could be used to put stringent constraints on the magnitude of these timing 
uncertainties. A thorough exploration of timing variation uncertainties for eclipse diagnos- 
tics is beyond the scope of the present analysis. Here, to evaluate simply the magnitude 
of these effects, we consider in some of our models the effects of arbitrarily time-shifting 
lightcurves produced by two different photospheric emission models, in such a way as to 
bring the mid-eclipse values closer to each other. This makes efforts to distinguish the mod- 
els more difficult, even though we are not formally fitting the eclipse center time in our 
analysis. In principle, the influence of a perturbing planet results in timing offsets which 
vary in amplitude from orbit to orbit. Here we are only considering the simple consequ ences 
of a constant, detrimental temporal shift in our models. Although lAgol et al.l (120051 ) have 
shown that a wide range of timing variations are possible in general (from seconds to minutes 
and more), we will focus here on a representative value of 7 seconds for definiteness, on the 
assumption that careful transit monitoring could plausibly constrain secondary eclipse tim- 
ing uncertainties to ^ 10 seconds in any given sy stem (which is about an order of magnitude 
improvement over current timing accuracies; e.g. iDeming et al.ll2006l ). 



3. The Method 



It is clear that current Spitzer data for a single eclipse in any of the three systems stud- 
ied here do not contain enough infor mation to determine the shape of the ingress or egress 
lightcurves with much accuracy (e.g. IDeming et al.ll2006l ). Multiple Spitzer eclipse observa- 
tions would be required to derive meaningful constraints on the day-side photospheric emis- 
sion properties of these hot Jupiters. The ability to distinguish between various photospheric 
emission models must then depend on whether or not the number of repeated observations 
needed is prohibitively large. To evaluate this, we simulate eclipse light curves produced 
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by different photospheric emission models, add noise in proportion to known measurement 
errors and statistically determine the number of observations needed to differentiate between 
two specific models. In other words, our analysis provides a model-dependent estimate of the 
level of atmospheric details that can be probed with repeated ingress/egress measurements. 



3.1. Simulating Eclipses 

To simulate eclipse light curves, a photospheric temperature model for the planet is 
orthographically projected onto a disk, which is then eclipsed using the known geometry 
of the system. At each point in time during ingress or egress, the emission from non- 
eclipsed regions on the planetary disk is integrated, in proportion to the (apparent) area of 
each visible disk surface element. Two analyzes are performed simultaneously. One focuses 
simply on the bolometric emission under the assumption that each surface element emits like 
a blackbody at the corresponding photospheric temperature. The other analysis is spectral: 
the contribution to the monochromatic flux at Earth is calculated for each apparent surface 
element by using the local photospheric temperature to interpolate in the above-mentioned 
grid of plane-parallel spectral emission models. Spectral contributions are then integrated 
over specific Spitzer response curves, or the wavelength ranges of the future JWST NIRSpec 
gratings, to determine the amount of light that would be collected by these instruments. 
The spectral models used are subject to strong systematic e ffects in specific spectral regions 



e.g. flux peak at ~ 4/zm, as seen Fig 1 of lSeager et al.ll2005l ). Testing how strongly our main 



results depend on such systematic spectral features is the reason why we use an alternative, 
very idealized bolometric blackbody emission model in some of our eclipse simulations. 

A 2D polar grid is constructed to describe the projected planetary disk, with 1000 radii 
sampled between and R p (the planetary radius) and 1500 polar angles sampled between 
and 27r. The radial resolution also determines the temporal resolution of our simulated 
eclipse curves since the stellar limb moves across the disk at a constant horizontal rate. We 
account for the curvature of the stellar limb, assuming it to be a perfectly opaque disk. We 
have tested our eclipse sim u lation apparatus by comparing results to the analytical solution 



described in I Warner et al.l (119711 ). With a [r,6] resolution of [1000,1500], a comparison to 
the analytical lightcurve solution for a straight stellar limb gives agreement at better than 
the 1% level (with growing errors towards the very end/beginning of ingress/egress because 
of grid under- sampling). The same level of accuracy is achieved in the case of a spherically 
curved stellar limb. 



In our analysis, we account for the influence of each system's specific geometry (i.e. 
inclination, ratio of planet and stellar radii, semi-major axis) on details of how the planetary 
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disk is being eclipsed. Table 1 lists the geometrical parameters adopted for the three systems 
of interest. The orbital velocity is calculated as v — 2ira/P, where a is the planet's orbital 
semi-major axis and P is the corresponding orbital period. The vertical offset (from the ob- 
server's point of view) between the centers of the planetary and stellar disks, h, is calculated 
as h = asin(90° — i), where % is the orbital inclination. The duration of ingress/egress is 
Egress = (\/(R* + R P ) 2 — h 2 — a/ (R* — R p ) 2 — h 2 )/v, where R* is the stellar radius. Figured] 
illustrates the eclipsing geometries for the TrES-1 and HD189733 systems. For simplicity, 
we ignore the effects of the small planetary rotational offset occurring between ingress and 
egress and consider instead that the sub-stellar point remains directly aligned towards the 
observer in all the models. The actual shift between ingress and egress is ~ 10° for the three 
systems of interest. We have performed explicit tests and found that such a small rotational 
offset of the photospheric emission model generally has little effect on our main results, as 
compared to other more important sources of uncertainty. 

As an example, Figure [2] shows the bolometric blackbody ingress and egress light curves 
obtained in the case of HD189733b for each of the photospheric emission models considered in 
our analysis. Dashed lines highlight the effects that strong planetary limb darkening has on 
the RadEq and Uniform models. It is obvious from these curves that the differences between 
the various photospheric emission models under consideration are subtle. It is also important 
to realize the existence of some level of degeneracy in our models. For example, during 
egress (only), there is little difference between the curves predicted by the Cho03-like and 
the Uniform models because the stellar limb exposes the Cho03-like photospheric emission 
features in such a way that contributions from hot and cold regions more or less cancel out. 
Finally, let us emphasize that the various photospheric emission models introduced above 
only define the shapes of ingress/egress lightcurves, in units of relative flux, in our modeling 
strategy. An absolute flux unit common to all these models is effectively selected only when a 
noise magnitude is chosen to generate a system- and instrument-specific stream of simulated 
data, as we shall now describe. 



3.2. Statistical Analysis 

With simulated ingress/egress light curves in hand, we can use known or estimated 
measurement errors to simulate noisy eclipse data. From these simulations, we can statisti- 
cally determine how much data is needed to differentiate between two specific photospheric 
emission models. For ease of comparison and definiteness, we chose the model without any 
heat redistribution, that is the RadEq model, to be our reference model. We generate an in- 
crementally increasing amount of simulated eclipse data for this model and determine when 
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it is that model lightcurves produced for other photospheric emission models become incon- 
sistent with the accumulated stream of data ( 1.6. clS db function of the number of repeated 
eclipse measurements). Our choice of the RadEq model as a reference is arbitrary, in the 
sense that any other model could be used as a reference for our statistical analysis. Notice, 
however, that by ruling out another model against the RadEq one, one effectively rules out 
some level of heat redistribution in the atmosphere studied. While any combination of two 
photospheric emission models could in principle be compared against each other, the interest 
of doing so is limited by the fact that all comparisons remain necessarily model-dependent. 
In the end, with our analysis we are only trying to assess the level of photospheric detail that 
could be distinguished by actual measurements. We will later suggest that non-parametric 
methods be considered for the study of hot Jupiter eclipses in the future. 

Throughout this analysis we work in units of relative flux. The planetary flux is nor- 
malized to its full day-side value, i.e. the value right before entering secondary eclipse 
(neglecting small rotational effects). A simulated data stream is generated on the basis of 
a cadence of observation and a Gaussian noise single-measurement error of variance o\ (in 
units of relative flux). The values of <7i adopted, as well as the cadence used to calculate 
the number of measurements occurring during one ingress or egress period, are estimated 
from previous measurements or simple extrapolations, as described in Section [3731 for Spitzer, 
and Section [5711 for JWST. Figure [3] shows an example of such a simulated data stream for 
HD189733b in one of Spitzer's IRAC bands. Noiseless ingress/egress curves for two pho- 
tospheric emission models (including the RadEq one used to generate the simulated noisy 
data) are over-plotted for comparison. Again, it is clear that repeated eclipses will be needed 
in this case if one wishes to discriminate between the two models, given with the level of 
noise in the data. 

The quality of any photospheric emission model in fitting a set of simulated data is 
evaluated by calculating the reduced \ 2 value for that model compared to the data. Since we 
are creating the simulated data from the reference RadEq model, with a sufficient number of 
repeated eclipse measurements, any other emission model will eventually provide a poor fit to 
the data (as measured by the reduced \ 2 value). The process of simulating data is repeated, 
as needed, for an increasing number of eclipses, N, with identical error properties but a 
random temporal offset added to each new eclipse set. This offset is added to account for 
the likely lack of accurate phasing even between successive eclipses in a real data set. For each 
value N of the number of accumulated eclipses, reduced \ 2 values are calculated comparing 
any photospheric emission model of interest to the data until, after enough measurements, 
models (different from the RadEq one) are ruled out at the 99% confidence level by the 
accumulated data set. 
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Figure H] illustrates this procedure. At each N, all models of interest are compared to 
the same set of simulated data, so that the random fluctuations in the data set affect all of 
the reduced x 2 values in a similar way, as is apparent in the fluctuation pattern common 
to each reduced \ 2 curve. The value of iV at which a reduced \ 2 nne crosses the 99% 
confidence level and does not return below it, iVo, is our measure of the number of eclipses 
needed to reliably differentiate between that model and the reference RadEq model[j] For 
example, Figure H] shows that the RadEq model can be differentiated from the CS05-like 
model with approximately 15 eclipses, since the latter is then ruled out according to our 
criterion. When a detrimental 7-second timing offset is added in the CS05-like model (in 
such a way that it makes harder to differentiate against the RadEq model), a total of ~30 
eclipses is required according to our criterion. A close inspection of Figure 0] suggests that 
this substantial difference could simply be a statistical fluke resulting from the use of one 
statistical realization of the data stream. 

To account for statistical variations in the value of Nq derived from single realizations 
of the data stream, we repeat the above procedure 500 times (i.e. we generate 500 data 
streams with identical parameters but different random number generators). From these 500 
realizations, we obtain a representative distribution of Nq values for a fixed set of model 
parameters. Figure [5] shows examples of such distributions of iVo values for Spitzer mea- 
surements. From these distributions, we now see clearly that a 7-second detrimental offset 
in timing applied to the CS05-like model does indeed make it harder to differentiate that 
model from the reference RadEq model. To be able to distinguish the fiducial RadEq model 
from the CS05-like model in 96% of the data stream realizations, 39 eclipses are needed. 
This number jumps to 52 eclipses if the detrimental 7-second timing offset is imposed on the 
CS05-like model curve, thus showing that eclipse timing uncertainties do affect the eclipse 
diagnostics being discussed here 

We note that other paramet r ic me thods, based on different eclipse diagnostics, ex- 



ist. In particular, IWilliams et al.l (120061 ) discuss the overall time shift between a Uniform 



m odel and an implementation of photospheric emission for the detailed atmospheric model 



of ICooper &: Showman! (120051 ). One advantage of our method over the time shift diagnostics 
is that we use more information, that is the full shape of the ingress/egress curves, which 
makes it less sensitive to degeneracies. To illustrate this point, we note that a comparable 
mid-eclipse time shift deviation from the prediction of the CS05-like model is expected for 
both the RadEq and the Uniform models, making it difficult to discriminate between these 



1 Alternatively, one could use the number of eclipses for which a reduced \ 2 l me first crosses the 99% 
confidence level to declare the corresponding model unfit. We have found that this less conservative criterion 
is unreliable. Figure [6] illustrates well how it fails. 
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two models solely on the basis of a time shift (see Fig. [2]). An adequately phased/rotated 
Cho03-like model could also induce a similar time shift relative to the RadEq and Uniform 
model predictions. Clearly, the issue of degeneracies is inherent to all parametric methods 
and we have already illustrated how our models also suffer from this limitation. In the future, 
non-parametric inversion techniques may thus be the best way to address this difficulty 



3.3. Estimating Spitzer Errors 



So far, we have considered in detail only the specific case of HD189733b observed by 
Spitzer at 8/xm with IRAC. The parameters used to produce the data stream in this case 
were inferred from actual observations of that system with Spitzer and from observations 
of another system with IRAC at 8/xm. Specifically, the single-measurement error was taken 
to be o"i = 0.31 and 134 measurements were assumed to be collected per ingress/egress 
period. Here we describe how we estimate these two numbers for every possible combination 
of instrument and planetary system of interest. This necessarily involves some level of 
extrapolation since many such observations have not been performed (or reported) to date. 

In the previous section, we defined <j\ as the error per observational point, for unity 
eclipse depth (i.e. in our relative flux units). In absolute units, this corresponds to the 
single- measurement error divided by the eclipse depth flux-equivalent, which can both easily 
be obtai ned from previ o us Sp itzer measurements when th ey exist. We have compil ed the 
results of beming et all (120051 ) using MIPS o n HD209458b ; Icharbonneau et all (120051 ) using 
the 4.5 and 8/im IRAC bands on TrES-1 and iDejninj^etalJ faoojl) using IRS on HD 189733b. 
The values of a\ derived from these measurements are shown in bold in Table 2. 

To estimate the errors for every possible combination of system and instrument, we pick 
an instrument, a specific spectral band, and rescale a known measurement error to derive 
the corresponding value for another system. Assuming that simple Gaussian noise statistics 
and the bright source limit are applicable, so that o> = yF, we expect the errors between 
system A and system B to scale approximately with the ratios of system fluxes as: 



&1,B 



*,B 



*,A 



p.A 
7 P,B 



where the * and p subscripts refer to the star and the planet, respectively. Using a simple 
blackbody-like scaling for the relative stellar and planetary fluxes, and isotropic emission in 
all cases, this expression becomes: 
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where d is the distance to the system, R is the object's radius, B\ is the Planck function, T 
is temperature, and A is taken to be the central wavelength of the instrumental band under 
consideration. The results of this extrapolation for o\ values are listed in Table 2. 

Table 2 also lists in each case the number of single data point measurements per 
ingress/egress period. These values are calculated from the duration of ingress/egress for 
each planet, combined with the cadence of each instrument. The ingress/egress durations, 
tgress, are listed in Table 1. The cadences for existing measurements are 12.3 s for MIPS, 13.2 
s for IRAC, and 14.7 s for IRS. New observing plans could result in different cadences, but 
these values should be sufficiently accurate at the level of detail of our analysis. We estimate 
that uncertainties in extrapolated o\ values may affect our results much more significantly. 



4. Results for Spitzer 

We performed an extensive survey of the parameter space of our eclipse models for each 
planetary system and Spitzer instrument of interest. The results of this exploration are 
summarized in Table 3. 

Figure M shows a typical reduced x 2 test for what we estimate to be the best scenario 
for Spitzer, differentiating between the RadEq and Uniform models at 8/im with IRAC in 
the HD189733 system. Clearly, Spitzer is not sensitive enough to differentiate between these 
two models in any practical te rmS | 

It is also impossible for Spitzer to differentiate the Cho03-like model from the RadEq 
model. One additional difficulty with the Cho03-like model is that the rotation of the 
temperature pattern would cause the detailed shape of the ingress/egress curves to change 
significantly over time if too many successive eclipses are required to reach good photometric 
accuracy. An instrument sensitive enough to differentiate the Cho03-like model from the 
RadEq model would probably have to do so in no more than a few successive eclipses to 
avoid excessive variations of ingress/egress shapes. In the specific Cho03-like model used 
here, the circumpolar vortices are ~10% cooler than the average photospheric temperature. 
At the level of detail probed by S pitzer, we find th at this model is indistinguishable from 



the Uniform model. (However, see ICho et al.l 120061 . for more extreme versions of this class 
of models.) 



5 According to our study, the only possible way to reliably discriminate between the RadEq model and 
the Uniform model with Spitzer, in less than a dozen eclipses or so, would be to discover a new system at 
least 10 times brighter than HD189733. 
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Of all the models we have considered in our survey, the only one that Spitzer could 
possibly distinguish from the RadEq reference model is the CS05-like model, in the bright 
HD 189733 system. Even in this case, from 9 to 52 eclipses are required (see Table 3 and 
Figure 4). This would be a very expensive observing program considering the strongly 
model-dependent nature of the diagnostics involved. 

If, instead of the spectral predictions of plane-parallel atmospheric calculations, we use 
simple, bolometric, blackbody photospheric emission models, we find in general that fewer 
eclipses are necessary to permit discrimination between two models (see comparisons in 
Table 3). Note that this difference is not due to the more efficient character of blackbody 
emission since all our ingress/egress models are eclipse-depth normalized in such a way 
that it is the normalized ingress/egress shape that is being used as a diagnostic. Rather, 
this shows that details of how the spectral energy distribution varies with local photospheric 
temperature do affect the diagnostics of our analysis. As additional observational constraints 
on hot Jupiter atmospheres become available, some classes of emission models (e.g. with or 
without clouds) will be favored over others. Nevertheless, the sensitivity of eclipse diagnostics 
to details of the spectral emission model assumed could remain significant. This is another 
important reason to favor diagnostics based on non-parametric inversion techniques in the 
future. 

While limb darkening could worsen the Nq statistics discussed so far, by reducing the 
information content of photospheric emission regions close to the planetary limb, we have 
not systematically explored this possibility in the context of Spitzer observations given the 
already poor statistics. We come back to this issue in the context of JWST observations 
below. 



5. The Case For JWST 

Since Spitzer will likely not be able to constrain photospheric temperatures and at- 
mospheric models at a very meaningful level with eclipse mapping, we turn our attention 
to the upcoming James Webb Space Telescope (JWST). The much increased sensitivity of 
JWST results from a combination of larger collecting area and reduced background. Of the 
instruments currently planned for JWST, we chose to focus on NIRSpec for our comparative 
analysis. This instrument will have three spectral gratings, covering the range of wavelengths 
from 1-5/xm. Spectra are not required for our analysis (or for eclipse mapping in general), 
but the dispersion of light provided by the gratings turns out to be useful in avoiding possible 
instrument saturation for the three bright systems of interest here. It may be possible to 
use the NIRCam instrument to study eclipses in dimmer systems than the three bright ones 
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considered here, thus increasing the total number of potential target systems, but we have 
not explored this possibility in detail. Our goal here is a first, very preliminary assessment 
of JWST capabilities in terms of eclipse mapping, relative to the Spitzer case. We caution 
that it is difficult to estimate high-contrast measurement errors for an instrument that was 
not designed for that task. 



5.1. Estimating JWST Measurement Errors 



As the JWST design is not yet finalized, our estimate of measurement errors for NIR- 
Spec is only approximate. We treat each NIRSpec spectral grating as a single bandpass, 
using the spectrally integ rated light only. We base our estimates on the detailed analy- 
sis of IValenti et al.l (120051 ) for a NIRSpec simulated (spectral) measurement of a secondary 
eclipse of HD209458b. The three gratings G140H, G235H, and G395H cover the wavelength 
ranges 1-1.8/im, 1.7-3/im, and 2.9-5/im, respectively. For HD209458b, the respective expo- 
sure times are estimated to be 2.4, 2.4, and 3.6 s. Adding a conservative estimate of 0.9 s for 
reset betwee n exposures results in cadences of 3.3 and 4.5 s for G140H/G235H and G395H, 
respectively (IValenti et al.ll2005l ). 



A calculation of the single-measurement error, entering our ingress/egress simula- 
tions is more complex and requires estimates of the number of stella r and planetary pho- 
tons that would be detected by each grating during one exposure. IValenti et al.l (120051 ) 
simulate noisy spectra for the HD209458 system, for two hours of exposure time. To es- 
timate the eclipse depth in any given grating, we divide the planetary spectrum by the 
stellar spectrum and sum over the wavelength range of that grating. Since the spectra are 
summed over two hours of exposure time, the number of photons collected during one ex- 
posure is easily obtained from the number of single exposures in a two-hour period. We 
assume that the nois e scales simply with the inverse square root of the number of detected 
stellar photons (like IValenti et al.l 120051 ). The single- measurement error is then given by 
(Ji = (a/A"* x N cxp )/(Np), where is the number of stellar photons, N p is the number of 
planetary photons a nd iV ex p is the number of unit exposures. Using Figures 3 and 4 from 
Valenti et al.l (120051 ) . a spectral resolution R = 2700 and a value of iV ex p based on the ca- 
dences listed above, we deduce the approximate single-measurement errors for HD209458b 
listed in italics in Table 2. Using the same system-to-system scaling method as before (§ 13. 3ft . 
we then extrapolate single-measurement <j\ errors for HD189733b and TrES-1 (see Table 2). 



As a consistency check on the <j\ values thus derived from IValenti et al.l (120051 ). the 
following simple arguments can also be used to scale values from the 4.5/xm Spitzer-IRAC 
band to the 4//m JWST spectral gratings. Wavelength ranges for these two bands overlap, 
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but they are not equivalent: the NIRSpec grating covers roughly 3-5/im, while the IRAC band 
covers approximately 4-5/zm. The additional 3-4/zm sub-range will allow JWST gratings to 
collect more photons and probe a region where the planet to st ar flux ratio is some what higher 
(due to the presence of a flux peak at ^ 4/zm; see Fig. 1 of ISeager et al.ll2005l ). Assuming 
that NIRSpec has 1000-fold increased sensitivity relative to IRAC, that the NIRSpec cadence 
is about three times that of IRAC, and that the JWST to Spitzer eclipse depth ratio in these 
bands is about a factor of two (from the 4/xm flux peak), we estimate a Spitzer-to-JWST 
or i ratio of ~50. This is consistent with the ratio 58 listed in Table 2. While this crude 
comparison lends support to the o\ values adopted below in our JWST analysis, there is 
clearly a need for better evaluations of the photometric accuracy expected with JWST for 
hot Jupiter eclipse mapping programs. 



5.2. Results For JWST 

The results of our survey of JWST/NIRSpec capabilities are summarized in Table 3. 
The much increased sensitivity allows relatively small photospheric emission features to 
be distinguished in the shapes of ingress/egress lightcurves in just one or a few eclipses. 
For HD189733b and HD209458b, the RadEq model would be easily distinguished from the 
Uniform and Cho03-like models by using the 4/xm grating. HD189733b is bright enough that 
any of the gratings would in fact be very useful, while for HD209458b, the 4/mi grating is a 
strongly preferable choice. The errors for this grating are expected to be significantly lower 
than the errors for the other two, partly because the planet to star flux ratio is higher at 
these wavelengths (according to the similar spectral models we and Valenti et al. 2005 used). 
A similar attempt to discriminate between the RadEq model and the Uniform model may 
require a prohibitively large number of eclipses for TrES-1. With the same 4/xm grating, our 
results indicate that the CS05-like model can be differentiated from the RadEq one on all 
three planets. We also find this time that, given the much increased sensitivity of JWST, 
applying a detrimental 7 s timing offset to the CS05-like model has little consequence for 
attempts to differentiate that model from the RadEq one. 

When we superimpose a strong limb darkening law on our Uniform photospheric emis- 
sion model, it can be more difficult to differentiate that model from the RadEq model. As an 
example, observing HD 189733b with the 4/zm grating, one eclipse is still enough to differen- 
tiate between the RadEq model and the limb-darkened Uniform model, but five eclipses are 
required to do so with the 1.4/im grating. At the level of precision achievable with JWST, 
it will thus be important to self-consistently account for limb darkening effects for reliable 
eclipse diagnostics. This is yet one more reason to focus on non-parametric inversion tech- 
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niques which, by construction, automatically incorporate any photospheric limb darkening 
present. 

As in Section HJ we find that the simplified bolometric blackbody analysis generally 
leads to a reduction in the number of eclipses needed for model differentiation, relative to 
the analysis based on the detailed spectral emission models. We interpret this as indicating 
that emission in the spectral models vary generally less with temperature (in relevant narrow 
spectral bands) than the Planck function does bolometrically. We note that a counterexam- 
ple of this general trend is found at the shortest wavelengths considered (the 1.4/xm grating). 
For this grating, the range of photospheric temperatures on HD 189733b conspires to make 
the emission properties in the spectral model vary more over ingress/egress than in the cor- 
responding blackbody model. This time, a lower number of eclipses is needed to differentiate 
between the RadEq model and the Uniform model, according to spectral models. Altogether, 
comparisons between spectral and bolometric blackbody versions of our eclipse analysis, and 
their qualitative agreement, suggest that our main conclusions (in terms of the number of 
eclipses required for model differentiations) are not too strongly affected by systematic ef- 
fects or uncertainties in the specific cloudless spectral emission models adopted throughout 
our analysis. 

As a special case, we have also considered the possibility of differentiating between the 
Cho03-like and the Uniform photospheric emission models on HD 189733b. We find that the 
4/xm NIRSpec grating should perform measurements precise enough to make the distinction 
between these two photospheric emission models possible with a single eclipse measurement 
(see the last three lines in Table 3). This would alleviate complications arising from temporal 
variations of the rotating atmospheric pattern, which are expected in Cho03-like atmospheric 
models. Let us recall that we have fixed the phase of the weather pattern in all our Cho03- 
like models to produce the greatest difference in ingress/egress shapes with respect to the 
RadEq model predictions. We selected a phase for the cold circumpolar vortices that would 
most easily be distinguished from that RadEq model. In a realistic situation, this specific 
phase would only be observed by chance coincidence. Our goal here was merely to determine 
the level of photospheric details that can be probed with JWST. This exercise shows that 
features of amplitudes and sizes similar to those of the circumpolar spots present in the 
Cho03-like model considered here (approximately 10% cooler than the average planetary 
photospheric temperature) are within the reach of JWST. 
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6. Summary and Conclusions 

The shape of ingress and egress curves produced when a hot Jupiter is gradually eclipsed 
by its parent star contains detailed information on the planet's photospheric emission prop- 
erties. Hot Jupiter atmospheric modeling efforts would greatly benefit from this information 
if it were to become available through applications of eclipse mapping techniques. In an 
attempt to evaluate the feasibility of eclipse mapping for hot Jupiters, we have generated 
ingress/egress lightcurves expected for a variety of plausible photospheric emission models 
and asked whether these models could be differentiated from one another, on the basis of re- 
peated eclipse measurements with current or future infrared facilities. We found that Spitzer 
has insufficient photometric accuracy to make efficient use of these eclipse diagnostics. Sim- 
ple scalings indicate that a system at least ten times brighter than HD189733b, currently 
the brightest known eclipsing hot Jupiter, would be needed for Spitzer to yield meaningful 
results in this respect. On the other hand, our preliminary analysis suggests that JWST 
has the potential to reveal relatively subtle features at the photospheres of currently known 
eclipsing hot Jupiters. 

Our analysis is voluntarily simple and limited in scope. By construction, we only com- 
pare specific photospheric emission models. This necessarily results in conclusions which 
are rather model-dependent. In addition, our treatment of atmospheric radiative trans- 
fer, and the related planet's photospheric emission, was relatively simple. We assumed the 
existence of a well-defined photospheric surface in all our spectral models, we combined 
simple horizontal photospheric temperature maps with the results of ID vertical radiative 
transfer calculations obtained independently and thus largely ignored important two- or 
three-dimensional geometrical effects in our analysis. Nevertheless, comparisons with crude 
bolometric blackbody emission models suggest that our main conclusions on the feasibility of 
using eclipse diagnostics with Spitzer and JWST are not too strongly affected by our various 
model specifics or shortcomings, because they mostly depend on the photometric accuracy 
achievable for the ingress/egress portions of infrared secondary eclipses. 

Throughout our analysis, we stressed the advantages of non-parametric inversion meth- 
ods for eclipse mapping, which would offer model-independent constraints on the photo- 
spheric emission properties of eclipsing hot Jupiters. By combining data at different wave- 
lengths, and thus probing photospheres at different heights, 3-dimensional emission maps 
of hot Jupiter atmospheres could potentially be obtained. It remains to be seen with what 
precision this concept can be applied in practice. It is likely that our model-dependent anal- 
ysis overestimates the level of details that can be probed by JWST with non-parametric 
eclipse mapping methods, but the promising nature of results obtained so far and the large 
potential value for the exoplanet science community provide strong motivations for a more 
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careful assessment of this exciting new possibility. 
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Table 1. System parameters adopted 





HD189733 


HD209458 


TrES-1 


distance (pc) 


19 a 


47 a 


150 b 


R* (R©) 


0.76 c 


1.146 d 


0.83 c 


T* (K) 


5050 c 


6000 f 


5250 b 


Rp (Rjnp) 


1.26 c 


1.35 d 


1.04 e 


T P (K) 


1200 


1440 


1160 


semi-major axis, a (AU) 


0.0313 c 


0.046 f 


0.0394 6 


orbital period, P (days) 


2.219 c 


3.524 f 


3.030 b 


orbital velocity, v (km s _1 ) 


154 


142 


142 


orbital inclination, % (° ) 


85.3 C 


87.1 s 


89.5 e 


star-planet offset, h (R p ) 


4.28 


3.62 


0.69 


ingress/egress t grcss (s) 


1770 


1520 


1050 



References. - iPerrvman et al.l ( 1997 ). nAlonso et al.l ( 2004 ). 
' Bonchv et a"D (l2005h . ferown et al.l (12001 ) fcozzetti et al.l (Eoo3), 
jMazeh et alLboOoh . fcharbonneau et all feoooh 
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Table 2. JWST and Spitzer single-measurement errors, <7i, and numbers of measurements, 

Ni, during a full ingress/egress 



Instrument, 


HD189733b 


HD209458b 


TrES-1 


wavelength (/xm) 






(71, iVi 


(7i, iVx 


JWST, G140H, 1.4 


1.5, 


536 


1.^, 461 


26, 318 


JWST, G235H, 2.3 


0.3, 


536 


0.^, 461 


4, 318 


JWST, G395H, 4 


0.02 


393 


0.ft£ 338 


0.3, 233 


Spitzer, IRAC, 4.5 


0.29 


134 


0.64, 115 


4.1, 80 


Spitzer, IRAC, 8 


0.27 


134 


0.71, 115 


3.8, 80 


Spitzer, IRS, 16 


0.46 


, 120 


1.3, 103 


6.2, 71 


Spitzer, MIPS, 24 


1.2, 


144 


3.5, 124 


17, 85 



Note. — Values in bold are calculated from previous mea- 
surement s and values i n itali cs are calculated from the esti- 
mates of IValenti et al.l (120051 ) . All other values were scaled 
accordingly (as described in Section I3TB1 . 



Table 3. Number of eclipses needed to differentiate between emission models with Spitzer or JWST. 



System 


Instrument, wavelength (/xm) 


Models compared 


N eclipses 


Spitzer 


any system 


any wavelength 


RadEq, Uniform 


>500 


HD189733 


IRAC, 8 


RadEq, CS05 w/no timing uncert. 


39 (78) 


HD189733 


bolometric w/8/xm noise stats. a 


RadEq, CS05 w/no timing uncert. 


9 (18) 


HDf89733 


IRAC, 8 


RadEq, CS05 w/7s timing uncert. 


52 (104) 


HD189733 


bolometric w/8/xm noise stats. a 


RadEq, CS05 w/7s timing uncert. 


10 (21) 


JWST 


TrES-f 


NIRSpec, any grating 


RadEq, Uniform 


>500 


HD209458 


NIRSpec, 1.4 or 2.3 


RadEq, Uniform 


>500 


HD209458 


NIRSpec, 4 


RadEq, Uniform 


2 


HD209458 


bolometric w/4/xm noise stats. a 


RadEq, Uniform 


1 


HD209458 


NIRSpec, 1.4 


RadEq, CS05 w/7s or w/no timing uncert. 


>450 


HD209458 


NIRSpec, 2.3 


RadEq, CS05 w/7s or w/no timing uncert. 


3 


HD209458 


NIRSpec, 4 


RadEq, CS05 w/7s or w/no timing uncert. 


1 


HD189733 


NIRSpec, 1.4 


RadEq, any model 


>500 


HDf89733 


NIRSpec, 2.3 


RadEq, Uniform 


>300 


HDf89733 


NIRSpec, 2.3 


RadEq, CS05 w/7s or w/no timing uncert. 


2(5) 


HDf89733 


bolometric w/2.3/xm noise stats a . 


RadEq, CS05 


3(11) 


HDf89733 


NIRSpec, 4 


RadEq, any model 


1(1) 


HDf89733 


NIRSpec, 4 


RadEq, Uniform w/ld b 


1(1) 


HDf89733 


NIRSpec, 4 


Uniform, Cho03 


2(2) 



Table 3 — Continued 



System Instrument, wavelength (//m) Models compared N ecKpses c 

HD189733 bolometric w/4/im noise stats. a RadEq, any model 1 (1) 
HD189733 bolometric w/4/im noise stats. a Uniform, Cho03 1 (1) 



a In the highly idealized "bolometric" models, ingress/egress curves are cal- 
culated assuming spectrally- integrated blackbody emission from the planet's 
various photospheric surface elements. Cadences and noise properties used 
to simulate eclipse data from these bolometric models are the same as the 
instrument-specific values assumed in all other models. 

b Model with additional limb darkening imposed. 

c The number of eclipses in parenthesis is for the value of HD189733b's radius 
recently revised down by 10%. See Appendix A for details. 
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Revised HD189733b Radius 



Recently, iBakos et al.l (120061 ) and lWinn et al.l (120071 ) have revised the value of the radius 
of HD189733b downward from 1.26 to 1.15 Jupiter radii. Accordingly, we have reconsidered 
our eclipse lightcurve analysis for this planet. The revised number of eclipses needed to dif- 
ferentiate between various photospheric emission models appears in parenthesis in Table 3. 
The smaller radius results in increased photometric errors and thus an increased number of 
required eclipses. Nonetheless, we find that the photometric accuracy of JWST at 4/im re- 
mains good enough for one or two eclipses to distinguish reliably among the various emission 
models considered. 
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Fig. 1. — Partially eclipsed temperature maps in a CS05-like model of TrES-1 (left) and a 
Cho03-like model of HD 189733b (right). The color-temperature scale (in K) is shown on 
the left of each panel. Notice how the different system geometries affect the orientation and 
shape of the eclipsing stellar limb, and consequently the detailed shape of the ingress/egress 
curves. 



-28- 




Fig. 2. — Simulated ingress and egress curves for different photospheric emission models 
of HDf 89733b: (a) RadEq model (black), Uniform model (dotted blue), CS05-like model 
(dashed red), and Cho03-like model (dot- dashed purple); (b) RadEq model (black), Uniform 
model (blue), RadEq model with strong limb darkening (dashed black), and Uniform model 
with strong limb darkening (dashed blue). Our analysis considers only the ingress and egress 
phases (separated by the vertical line) and thus excludes the full occultation period. 
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Fig. 3. — Example of ingress/egress simulated data for HD189733b in Spitzer's 8/zm IRAC 
band, for single-measurement noise of variance G\ = 0.31 (crosses) based on existing Spitzer 
measurements. For comparison, two noise-free curves are shown for the RadEq model (black) 
and CS05-like model (red dashed). Here again, our analysis excludes the full occultation 
period. 
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N (number of eclipses) 

Fig. 4. — Reduced x 2 values of three different model fits to a simulated dataset for 
HDf 89733b, as a function of the number of repeated eclipse measurements, N. The single- 
measurement noise variance adopted in the simulated ingress/egress dataset is <7i = 0.31. 
The 99% confidence limit above which models are ruled out is shown as a dashed line. The 
RadEq model used to generate the simulated ingress/egress dataset is shown as a black line. 
With a large enough set of eclipse data, other models (red line: CS05-like model; blue line: 
CS05-like model with detrimental 7 second timing offset) are eventually ruled out at high 
significance. 
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Fig. 5. — Histograms of iVo, the number of eclipses needed to differentiate at the 99% 
confidence level between fits to the correct underlying RadEq model and a CS05-like model 
(a), or a CS05-like model with an imposed 7 second detrimental timing offset (b). These 
distributions emerge from 500 independent realizations of repeated eclipse sequences like the 
one shown in Figure HI Nq is defined as the value of N in such sequences at which each 
model's reduced x 2 curve reliably crosses the 99% confidence level. 



-32 - 




Fig. 6. — Example of the reduced x 2 values as a function of the number of eclipses, N, for 
data compared to the RadEq model (black) and the Uniform model (red) for HD189773b 
using Spitzer's 8/im IRAC band. The models cannot be differentiated by this method. 



